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1. Introduction 



In order to compute the spectrum of a self-adjoint operator A in a Hilbert space H one can 
approximate A by a sequence of finite-dimensional operators An = P„AP„, n = 1, 2, . . . , where P„ 
I— I are orthogonal projections on finite-dimensional subspaces of H with the property P„ — >■ / in the 

strong sense as n ^ oo. However, it is well known that the operators An may have eigenvalues 
which in the limit do not converge to spectral points of A. Such eigenvalues are termed spurious 
eigenvalues and their presence is often described as spectral pollution (see e.g. [3], [5]). The 
"j^ spectral pollution usually happens in spectral gaps of the operator A. This phenomenon has been 

extensively studied both in the abstract setting (see e.g. [7], [1]), and in various special cases (see 
I I e.g. [2], [S], [6]). In particular, it was shown in [7], [4], [9], and [5] that the spectral pollution may 

occur at any point in a gap of the essential spectrum of A. 

Perhaps, the simplest example illustrating the spectral pollution, is the classical Toeplitz ma- 
trix. Let a be a real- valued piecewise continuous function on the interval [— tt, tt), and let P„ be the 
\^ orthogonal projection in i^(— 7r,7r) on the subspace spanned by the exponentials (27r)^^/^e^*'^^, 

fc = 1, 2, . . . , n. We define the Toeplitz operator with the symbol a as 

• Tn = Tn [d] = Pn APn ; 

where A is the operator of multiplication by a. If the range of the function a is a disconnected 
set, then the spurious eigenvalues, in the limit n — ?► oo, fill in the gaps separating the components 

^-H of the range. More precisely, it can be inferred from [1] that an open interval I strictly inside the 

^ gap contains Wlogn + 0{1) eigenvalues of T„, where W = W{a,I) is an explicitly computable 

• '-j constant. The objective of this note is to study in more detail the spectrum of the Toeplitz matrix 

rS for the piecewise constant symbol of the form 

, . , . ( 0, X e [-TT, L), 

^'■'^ "("^ = 1 1, ^e[L,n), 

where L e [0,7r). In !5j, it was shown numerically that if Lp ~ nq with some (co-prime) integer 
numbers p,q, then the spurious eigenvalues of the operator r„[a] are "nearly periodic" in n with 
a "period" w = uj{p,q), see formula (2.1) below. Specifically, let A^*^-* G (0, 1) be an eigenvalue of 



the matrix Tfc[a]. Then for each of the operators T/t+tj/ [a] , I = 1; 2, ... , there exists an eigenvalue 
such that the difference aC'+'^^'+i)) - A(''+'^') tends to zero as ? ^ oo. Ahhough the 
convergence rate of the latter was not estimated in [5], graphs of the eigenvalues qualitatively 
showed a rate which is faster than the logarithmic filling rate of the gap. As an example. Fig. [T] 
(see [5]) gives a diagram of the numerically computed eigenvalues of P„AP„ vs. n, with L = ^. 
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The symbols "o" and "x" mark the eigenvalues for even and odd n respectively. The diagram 
clearly suggests periodicity with period 4. 
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Figure 1. Eigenvalues of T„[a], with L = |, plotted vs. n. 



The paper ^ (see also llOj in this volume) offered no quantitative measurement of the properties 
of the periodic behavior. Moreover no rigorous proof has been given, see |10) for some intuitive 
discussion of the phenomenon. 

In this article we address the question of the near periodicity of eigenvalues, but instead of 
the symbol it is more convenient to take a of the form (2.1). We have not been able to find 



a proof of this effect for the operator T„[a] itself, but we c an s how its presence for the squared 



Toeplitz operator, i.e. for M„[a] = (T„[a])^, see Theorem 2.1 Numerical examples and more 



detailed conjectures with regard to the periodicity will be presented in a further publication. 

2. The main result 

We are concerned with the spectrum of the squared Toeplitz operator 

M„ = A/„[a] = (T„[a])2. 

First, we introduce some consistent notation for eigenvalues of various operators which appear 
later in the paper. For any self-adjoint matrix operator 5* = S'„ of size n, we shall denote by 
Xj{S), j = 1, . . . , n, its eigenvalues labeled in the descending order, and by Hk{S), fc = 1, . . . , n, its 
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eigenvalues labeled in the ascending order, so that 

flj{S) = \n+l-j{S) . 

Also, for brevity, in the particular case of operators Mn, we shall write 

The next theorem is the main result of the paper: 
Theorem 2.1. Let a{x) be the function 



(2.1) a{x) 



-1, X e [-TT, L), 
1, X £ [L,tt). 



with some L g [0,7r). Suppose that 

(2.2) pL = nq 

for some co-prime p £ N and g G Z. Define 

'2, q = 0, 
oj = uj{p, q) = { p, p and q are odd, 

2p, either p oi q is even. 

Let e e (0, 1) be a fixed number, and let j be such that /uj"'' < 1 — e. Then for a sufficiently large 
number K > and n > Kuie~^ , 

(2.3) |M"^-M"+")|< ^"^^ + ^°g'"\ 

■' en 

with a constant C > independent of p, q,n, j,e. 

In general, by C and c we denote various positive constants independent of n,p, q and e, whose 
precise value is of no importance. 

Note that under the condition (2.2 1 we have 



(2.4) uj{tt + L) = mod 2tt 

mmediately implies that u 
and all sufficiently large n 



Remark 2.2. Theorem 
any fixed fixed N — 1,2, 



2.1 



(n) 

immediately implies that under the condition ^ - < 1 — e, e > 0, for 



^2 5) ^^{n+uj{m+l)) _ ^(n+L.m)| ^ ^^(1 + log^ n) 



en 

for all m = 0, 1, . . . ,iV, uniformly in j. This means that for any N the spectra of M„[a] have 
strings of length N of "nearly equal" eigenvalues. 



Throughout the proof of Theorem 2.1 wc omit the subscript n for brevity whenever possible. 
Denote Q = Qn = I — Pn and compute, remembering that = L: 

M„[a] = {PAPf ^ PAPAP = PA^P - PAQAP = P - PAQAP. 

This means that 

mJ"' = 1-A,(i?„), 



where = Bn[a] — PAQAP. Therefore it suffices to prove the inequality (2.3) for the eigenvalues 
Xj{Bn) which satisfy Xj{Bn) > e, instead. 
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The entries of the matrix i3„ are easy to find: 



(2.6) 



6w = ''w = 7— ar-mttm^i, r,l ^ 1,2, ... ,n, 



m<0, 
7n>n+l 



where the Fourier cocfiicicnts 
are given by 

(2.7) flfc = 



V27r J- 



2 (-1)^- 
TT ik 



k = 0, 



[1 - e*'=(^+'^)] , k^O. 



The crucial point of our argument is that due to (2.4) the exponential on the right-hand-side of 



(2.7) is w-periodic as a function of k. This fact is used only once, in the proof of Lemma 2.6 
We need to establish the following theorem. 



Theorem 2.3. Let a(x) be as defined in (2.1), and suppose that the condition (2.2) is satisfied. 
Let e G (0, 1) be a fixed number, and let j be such that \j{Bn) > e. Then for all sufficiently large 
n> I, 



(2.8) 
uniformly in j . 

First we estimate the entries fei?'' 



|Aj(-B„) — Xj{Bn+uj)\ < 



Ca;(l + log2 n) 



Lemma 2.4. Let a{x) be defined by (2.1), and let b^^l^ be defined by (2.6). Then for all \ < I < 
r < n we have: 



(2.9) 

and 

(2.10) 



(„) 16 1+logn 



Proof. Substituting ^J} in ( |2.6| we get 

E 



1 



(2-11) = ^ , w 

m<0, ^ ' 
m>n-t-l 



Therefore 



|6w|<-, E 



m>n+l 



1 



(r — m)(m — I) 



1 



4 

Tr^\\r-{n + l)\ |(n + l)-/| ' rl 



E 



m< 
m>n-|-2 



^ |r — m| |to — Z| 
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The second term can be easily estimated by an appropriate integral. If r > Z, then 



E 

■m< — 1, 
m>n+2 



1 



r — m m ■ 



< 



1 



(2.12) 



1 



(r — x){l — x) 
log 



dx 



1 



+1 {x~r){x-l) 



dx 



n + l-l 
n + 1 — r 



log 



The bound (|2.9|) follows immediately. 

]. Ii 

log 



Let us derive (2.101. In view of the straightforward bounds 

+ i-r 



n + 1 — r 



< 



-I 



n + 1 — r 



log(y) 



< 



- 1 



the right hand side of (2.12) does not exceed 

1 1 

n + 1 — r I 

For the case r = I, a similar estimate can be obtained: 



E ■ 

m < — 1 , 
m>n+2 



Hence for all r > I we have: 



1 



< 



1 



_i {x — r)^ 
1 1 



dx 



(r — x)' 



-dx 



n + 1 — r 



Ifowl < 



1 



TT^ \ n + 1 — r 



which coincides with (2.10). 



Let w e N, 1 < w < n, and let 



n+1 



n even, 



n odd. 



□ 



Let us construct, using matrix Bn+^i, two new auxiliary matrices. The (?i + cj) x (n + a;)-matrix 
Dn+co bas the entries 



(2.13) 



,(n+cj) 



k<r <k + u}, or k<l<k + ij, 







otherwise. 



The n X n matrix has the entries 



(2.14) 



, (n+w) 
, (ri+w) 
, (n+w) 



l<r<fc-l,l</<fc-l, 
k<r<n,k<l<n, 
\<r<k~\,k<l<n, 
k<r <nA<l<k-l. 



The method of constructing the matrices D^^^ and F„ from B^+uj is illustrated by Figure|2] First, 
we shade the central "cross" of uj rows and uj columns in B^+u), starting with row and column 
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number k. The matrix Dn+u is eonstructed by replacing everything outside the cross by zeros, 
and the matrix by "removing" the cross, and pulhng the remaining four blocks together. 



/ ^1,1 • 


b\.k-\ 


buk ■ 


b\,k+co-l 


b\,k+co 


b\.,i+Q} \ 


bk-\,\ 


bk-\.k-\ 


bk- 1 ,k 


bk-\.k+a-\ 


bk-l.k+co 


bk-l.n+m 




bk,k-l 


bk,k 


bk,k+CD-\ 


bk,k+a 


bk,n+(0 


bk+(o-\,\ ■ 


■ bk+CD-l.k-\ 


bk+co-\,k ■ 




bk+m-l.k+ai ' 


bk+0)-l,n+co 


bk+a,] 


bk+a.k-\ 


bk+0),k 


bk+co.k+(o-\ 


bk+a.k+co 


bk+m.n+d) 


\ bn+coA 


bn+co.k-l 


bn+a),k 


bn+co,k+a)—l 


bn+co.k+co 


bn+co.u+o) J 



Figure 2. Matrix 3^+^ and the central "cross" used in the construction of ma- 
trices Dn+uj and F„ 

Introduce the projections Un^uj = Pk+uj-i ~ Pk-i, and En,uj = Pn+uj{I - ^n,uj)- In terms of 
these projections, the operator can be represented as follows: 

P^n-\-u Iln.aj-^n+cjIIri,^ ^ ^n.uj -E^n-\-uj^n,uj ^ ^n.ijj Pn-\-uj^n,uj ■; 

SO that 

The next lemma estimates the matrix Dn+p in terms of these projections. 



Lemma 2.5. Let 6j," be defined as in (2.6|, and c?|," ^ as in (2.13). Then for any e > and all 

n > p we have 



(2.15) 



\{D,,+pU,u)\<ie + Cpn ^)||n„,pu|| 



Cp 1 + log n 2 



for all u £ L^(— TT, tt). 

Proof. Write the straightforward estimate: 

< ||ii„,^B„+„n„,^|| ||n„,^u|p + 2||n„,„B„+^s„,^M|| ||n„,„u|| 
(2.16) < (||n„,^B„+„n„,^|| + e)||n„,„u|p + e-i||n„^^s„+„s„,^||2||s„,^w|i2. 

Here we have used the elementary estimate 2ab < ea^ + e~^6^, e,a,b > 0. 

Let us estimate the matrix norms entering the above inequality. For the first term we estimate 



the Hilbert-Schmidt norm using the bound (|2.10) 
(2.17) 



16 



Iln.Lj-Bn+wIIn.LjIle, ^ 2- 



2 k+ui—l r 



E E 

= k l=k 



n + uj + 1 — r I 
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For the matrix ||n„,^_B„+;^^„_(j|| we also estimate its Hilbert-Schmidt norm, but now we need both 
(pjl) and (pl0|; 



(2.18) 



k-\-Lj—l k — 1 n-\-uj — 1 

r—k l — l r—k+uj l—k 



^(ri+w)|2 



Let us estimate the first sum, which we denote Si. Spht it into two parts: 1 < Z < fc/2 and 
k/2 < I < k — 1. For the first part we use the estimate (2.9), so that 

|^(„+.)|^^l + logn^ 



and 



k+uj-l 



/111 ^2 k/2 



(1 + logn)2 



r=k l<k/2 



r=k 1=1 



The part with k/2 < I < k — 1 is estimated with the help of (2.10), so that 

C 



,(n+cj)| 
^rl I 



< 



and 



Thus 



k+u)-l 

E E 1^: 

r=k k/2<l<k-l 



Si < CUJ 



(1 + logn) 



The same bound holds for the second sum on the right hand side of (2.18). Together with (2.171 
and (2.16) these bounds lead to the claimed estimate (2.15). □ 



Lemma 2.6. Suppose that (2.2) is satisfied. Then for all n > 1 we have 



(n) (n) 

Proof. First we estimate the difference /^^ Ki ■ convenience we re- write the formula 



(124 for 5'';+"^ 



- E 

27r ^ 



m<0, 
7n'>n-\-uj-\-l 



(2.20) 



rn<0 m>n+l 



Our plan is to estimate the Hilbert-Schmidt norm of - So, we estimate carefully the 

sum of squares \ fri — bri\'^ over the four ranges of r and I, specified in (2.14). 

Case 1: l<r<fc-l, l</<fc-l (upper-left block in Figure [2). According to ( |2.20 | 
we have 

f (") _ _ J_ \ " / ^ — n n \ 

J rl rl o / ^ V ' — m—uj^ra+u—l "r— m'^m — / j • 



m>n-t-l 
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Using ( |2.7[ ) and ( |2.4 ) we get the bound 



-1^ 



1 



(r — 771 — uj)(77i + Lu — I) (r — m)(m — /) 



< 



8uj 



1 



(r — m)2(TO — Z) (m — r)(m — Z)^ 



Therefore, we can estimate: 



rn>n+l 



1 



(r — m)2(TO — (m — 7')(m — Z)^ 



(n + 1 — rj(n + i — <j 



and hence 
(2.21) 



k-l 



i,r=l 



Case 2: k<r< 77, k<l<n (lower-right block in Figure [2]). Using again (2.20), from 



( |2.14[ ) we get: 

f (") _ _ J_ \ " / „ _ ^ /, 

m<0 

Arguing as in Case 1, we obtain the same bound: 
(2.22) 



in) .(«)|2 ^ CUJ^ 



Lr—k 



The remaining two cases are trickier: 

Case 3: l<r<k—l, k<l<7i (upper-right block in Figure [2]). By (2.141, 



I rl Ir r,l-\-uj Ir 

1 

2tj: 



27r 

m>n+l 

Arguing as in the previous case, we obtain 



•n<0 



l^r—m—Lu^m — l ^r—m^m—l\ — 



1 



TT (r — m)(m — Z)^ 

8w 1 
TT (r — m)2(m — Z) 



, m < 0, 
, m > n + 1. 



We estimate 



Sri 



E 

m<0 



1 



(r — TO)(m — Z)-^ 



< C 



1 + log 71 
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and 



Therefore 



E 



m>n+ 
fc— 1 n k—1 



^ (r — m)^{m — I) 



< c 



1 + log n 



r=l l=k r=l l=k 

Combining these bounds we arrive at the estimate 

k—l n 

(2.23) EEI/h^-^w^I'^ 

r=l /=fc 



C(l + log' n) 



in) ^(n)|2 ^ Ca;'(l + log n) 

n 2 



Case 4: fc<r<n, 1<^<A; — 1 (lower- left block in Figure [2]). Since the matrices F^") 
and are Hermitian, we can use the estimate obtained in Case 3, and hence 



(2.24) 



n k — 1 

EEi/^, 



(») ^(»)|2 < Ca;2(l + log'n) 



End of proof. We estimate the norm by the Hilbert-Schmidt norm: 

\\Fn — Bn\\ < \\Fn ~ i?„||62- 



In view of (2.21), (2.221,(2.23), (2.24|, 



\\Fn-Bn\\62 < Cu 

which coincides with the proclaimed estimate. 



1 + log n 



□ 



Our next step is to estimate the difference between eigenvalues of Bn+^ and i^„. Instead of 
the eigenvalues themselves, it is more convenient to work with their counting function. Denote by 
n+(A, S"), A > 0, the number of eigenvalues of a Hermitian matrix S strictly above A. 

Lemma 2.7. Suppose that A > and n > KcoX^^ with a sufficiently large constant K > 0. Then 



n+[\ + 



Co; 1 + log n 



Fn < n+{X, Bn+uj) < n+ A 



Coj1 + log n 



An) \ A n 

Proof. We use the projections S„.^,n„.^ introduced earlier. Recall that 



1 n 



Bn-\-uj — ^n,LjBYi-\-uj^n,uj ^" -^n+oj ■ 



By Lemma |2.5| this implies that 



Cwl + log'n„ ^ -iNT-r 
^n,uj + (d + Cujn )II„^^, 



Bn+Lj ^ ^n,LL>Bn+ui^n.ui f ^n,u} — ((5 + CiOn "'")n„.^, 

n 

for any S > 0. Thus the counting function n+(A,-B„+(^) satisfies the following estimates: 

Cuj 1 +log'n. 



(2.25) n+{X,B„+u) < n+[ Bn 



(2.26) 



n+(A, Bn+uj) ^ ^1 ^n.,uiBnJrUj^n,. 



^n,u + n+{\ {5 + Cujn )n„^„). 



Clo 1 + log^n„ 
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with an arbitrary S > 0. Now take 6 — A/2, so that for sufficiently large K under the condition 
n > KujX^^ we have 5 + Cijjn^^ < A. Therefore the second term on the right-hand-side of (2.251 
equals zero. The matrix Fn is obviously similar to S„^^i?„+tjS„^tj, so that their positive eigenvalues 
coincide. Thus (2.25) and (2.26) lead to the required inequalities. □ 



Proof of Theorem 12.31 By Lemma [2^ the elementary perturbation theory yields: 



for any A € (0, 1]. Using this bound in combination with Lemma 2.7 we get 

Cwl + logSi \ ^XR^^ A Cujl+}o^ 

AH ,Bn+uj < n+{X,Bn) <n+[\ , S„+„ 

en J \ e n 

for all A G [e, 1] and n > Koje^^, with a sufficiently large constant K. This means that if Aj(i?„) > 
e, then 

en 

for all n > i^we^i. □ 

As we have already pointed out. Theorem |2.1| follows from Theorem |2.3| due to the equality 
A^^"^ = 1-A,(S„). 
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